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Sublattice ordering in a dilute ensemble of defects in graphene 
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Abstract. - Defects in graphene, such as vacancies or adsorbents attaching themselves to carbons, 
may preferentially take positions on one of its two sublattices, thus breaking the global lattice 
symmetry. This leads to opening a gap in the electronic spectrum. We show that such a sublattice 
ordering may spontaneously occur in a dilute ensemble defects, due to the long-range interaction 
between them mediated by electrons. As a result sublattice-ordered domains may form, with 
electronic properties characteristic of a two-dimensional topological insulator. 
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The flatland of graphene [T] turned out to be a lab- 
oratory where fundamental ideas of particle physics and 
cosmology [2| find an on-table realization. The low-energy 
spectrum of charge carriers in graphene resembles that of 
massless relativistic particles [3]: it directly follows from 
the symmetry of its hexagonal lattice [3~J[4] that the spec- 
trum is gapless and linear. In pristine graphene, the lattice 
symmetry is so robust and resilient against spontaneous 
symmetry breaking that no gap-generation mechanism has 
yet been proposed. Here we show that a relatively low con- 
centration of adsorbents sprinkled over graphene flake can 
reverse the situation. It turns out that even at relatively 
low concentration the adsorbents may self-organize into 
the partially ordered state. The ordering generates a gap 
in the electronic spectrum and converts graphene into a 
topological insulator. 

The honeycomb lattice of graphene with two carbons in 
the unit cell can be considered as a combination of two 
overlaying triangular sublattices, A and B - see Fig. 1. 
It is the symmetry between the sublattices that prevents 
formation of the gap in pristine graphene [3j[4] . Consider 
now one-sided chemisorption of atoms on a flake deposited 
on an insulating substrate. Some adsorbents, such as flu- 
orine, hydrogen, or hydroxyl group are known to form a 
covalent bond with a particular carbon [H]-[8], thus violat- 
ing the symmetry between the sublattices locally. As long 
as the density p of the adsorbents is low (p <C a~ 2 , where 
a is the lattice constant of graphene) and their locations 
are random, one should expect approximately equal occu- 



pation of the sublattices, pa ~ Pb (here, Pa/b is density 
of adsorbents on A/B sublattice, and p = pa + Pb)- Ef- 
fect of the local symmetry violation is, then, limited to the 
scattering of electrons. 

However, under certain conditions the global sublat- 
tice symmetry can be spontaneously violated (pa 7^ Pb), 
leading to the correlations between sublattice affiliations 
of adatoms. These correlations between adsorbents' po- 
sitions are electron-mediated: each adsorbent creates a 
sublattice-asymmetric perturbation in the electronic den- 
sity (which is equivalent to the polarization of the electron 
states inside the unit cell), which in turn interacts with 
other adsorbents. The effective long-range pair interaction 
between the adsorbents depends on whether they reside on 
the same sublattice, or on different ones. This situation 
can be described in terms of the two-dimensional Ising 
model with the direction of an Ising 'spin', m = ±1 iden- 
tifying to which sublattice each given adsorbent belongs. 
The finite difference pa~Pb = pM between the sublattice 
occupations emerges, then, as an analog of the sponta- 
neous spin polarization M = (m) in the Ising model, and 
spontaneous breaking of sublattice symmetry over meso- 
scopic scale areas in a flake gives rise to the gap in the 
electronic spectrum, Fig. 1. 

Effect of the partial ordering illustrated in Fig. 1 is de- 
termined by the relation between the electronic de Broglie 
wavelength A = h/p and the mean distance, p -1 / 2 between 
the adsorbents. While 'high-energy' electrons with short 
wavelength A -C p^ 1 ^ 2 are insensitive to the ordering (for 
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them the adsorbents act as individual scatterers), for elec- 
trons with long wavelength A > p^ 1 / 2 the ordering implies 
global symmetry-breaking field, which opens the spectral 
gap, A [see Fig. 1]: 

e = ±y / v 2 p 2 + A 2 , A = Mua 2 p. (1) 

Here, v is Dirac velocity of electrons, and u is the elec- 
tron - adatom coupling constant. The density of states 
(DoS) corresponding to the gapped spectrum Eq. 
is 7 = 7o#(£ — A) (here 9(z) is the step function and 
7o = 2e/irv 2 h 2 is DoS in pristine graphene, which takes 
into account the spin and valley degeneracy). One can 
say that in a unit area of the graphene flake as many as 
A 2 /(ttv 2 H 2 ) states that originally resided within the gap 
vp < A have been transferred by the ordering to the higher 
energies, \e\ ~ hv^fp. In the absence of carriers the con- 
duction band is empty, while the valence band is full. The 
gap formation thus reduces the sheet density of the en- 
ergy by, approximately, A 2 yfp/hv ~ Al 2 u 2 (a 2 p) 2 (y/p/hv) 
For the ordering to be thermodynamically favorable, the 
energy gain should exceed the free energy cost, kgTM 2 p 
related to the entropy drop upon ordering. We therefore 
predict a phase transition with the sublattice symmetry 
breaking at T = T c , 

T c = C^{a 2 pfl 2 (2) 
hv/a 

where C ~ 1 is a numerical factor which will be evaluated 
below. The formation of the gap in the electronic spec- 
trum due to the partial ordering of adsorbents [see Fig. 1] 
should manifest itself in ARPES and in a suppressed in- 
frared light absorption. 

To describe the sublattice ordering transition we use the 
RKKY-like approach |9j to the pair correlations between 
adsorbents (adatoms) and map the problem onto the Ising 
model. The Hamiltonian of graphene with adsorbents re- 
siding on carbons has the form [8|[l0] 

H = va ■ p + ua 2 ^[a z m l + V Ies ]S(r - rj) (3) 

i 

The first term H describes free electrons with linear spec- 
trum. The second term accounts for the interaction of the 
electrons with the adsorbents. It is responsible for the A- 
B symmetry breaking. Three Pauli matrices, <J x ,y,z act on 
the sublattice components of the electronic Bloch function. 
Ising 'spin' m, determines the sublattice A (ra, = +1) or 
B (rrii — —1) that hosts given adsorbent. The residual 
term V^. os contains two contributions, which do not affect 
sublattice ordering. Firstly, there are two electronic val- 
leys in graphene, and V Ies takes into account inter-valley 
scattering. Although we dropped the valley indices in Eq. 
([3]), everywhere below the valley degeneracy is taken into 
account. Secondly, V rcs contains a "scalar" contribution 
describing the sublattice- and valley- independent channel 
of electron scattering. Although such a contribution does 




Fig. 1: (a) Disordered (top) and sublattice-ordered (bottom) 
state of a dilute ensemble of adsorbents residing over the sites 
of the honeycomb lattice of graphene. (b) The 1st Brillouin 
zone of graphene (dashed lines), with three reciprocal lattice 
vectors linking triads of equivalent corners, such as points K. 
The top/bottom insets illustrate that partial sublattice order- 
ing opens a gap in otherwise gapless electronic spectrum of 
graphene. 

not affect sublattice ordering directly, it is helpful to notice 
that it leads to a repulsive interaction between adatoms 
and precludes adatom clustering discussed in 

A remark is due here concerning the Hamiltonian |3]). 
Strictly speaking, the interaction term in the ^ is not 
well defined due to ultraviolet problems and is written 
in this form for illustrative purposes. In a more rigorous 
approach a z mi + V Tes should be replaced by an energy- 
dependent T-matrix defining the long-distance asymptotic 
form of the electron wave scattered off the defect at point 
Yi . The typical energy of electrons involved in the RKKY 
exchange between impurity atoms is Ae ~ hv^fp. If T(e) 
is a slowly varying function of energy for |e| < Ae, one 
can neglect its energy dependence altogether. Technically, 
this amounts to treating the scattering terms in Eq. Q 
in the Born approximation. Under special circumstances, 
for instance in the presence of resonant impurity levels at 
energies comparable to Ae, such an approximation breaks 
down and an alternative approach is to be used, e.g. such 
as in [11] . 

An adsorbent l V attached to one of carbons redis- 
tributes the electron density between A and B sites. The 
density redistribution creates a sublattice-asymmetric po- 
tential landscape for another adatom £ j' at a distance ry. 
We have found that energy of two adatoms is lower when 
they reside at the same sublattice than when sublattices 
are different: 

i¥=3 3 

and Tij = \vi — Yj\. The electron- mediated interaction 
E{ nt between adatoms can be derived using the standard 
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diagrammatic technique [9]. The interaction Eq. Q re- 
sembles the RKKY interaction between spins in dilute fer- 
romagnetic semiconductors [l2j. There is however an im- 
portant difference specific for undoped graphene [13] : cou- 
pling E- mt monotonically decreases as r~ 3 , while RKKY 
coupling oscillates and decays in two dimensions only as 

y 

Conventional RKKY interaction is caused by Friedel os- 
cillations - the electron density around a point defect oscil- 
lating in space with the wave vector equal to the momen- 
tum transfer upon scattering between the extreme points 
of the Fermi surface. In graphene the hexagonal Brillouin 
zone (BZ) is characterized by two triads of equivalent cor- 
ners, K and K' related within each triad by the reciprocal 
lattice vectors, G, G', G" [see Fig. 1(b)]. In the absence of 
charge carriers the Fermi line is shrank to the points K and 
K' which become its extremes. Perturbation of the elec- 
tron Vanier functions by an adsorbent creates oscillating 
charge redistribution analogous to the Friedel oscillations. 
The wave numbers in this case are equal to the reciprocal 
lattice vectors, e.g. G,G', or G". Usually the amplitude 
of Friedel oscillations is proportional to the Fermi level 
DoS of electrons. The latter, formally, vanishes: in neu- 
tral graphene 7 cx e = vp, e.g., 7(0) = 0. Nevertheless 
the charge density modulation still takes place. Indeed, 
according to the Heisenberg uncertainty principle, given 
the distance between the two adsorbents, the charac- 
teristic electron momentum counted from the BZ corner 
is p ~ h/rij. Therefore, for the estimation of RKKY cou- 
pling we should use 7 ~ (fiwy) . This is the source of 
the additional l/r^ factor in the RKKY interaction be- 
tween spinfull on-site impurities in graphene fl3| and in 



Eq. Q describing the distance dependence of the effective 
interaction between adsorbents. 

Equation Q determines the Ising model in the ensemble 
of randomly distributed 'spins' with a long-range interac- 
tion. To evaluate the critical temperature T c we simulated 
the transition numerically: we used 10 realizations of ran- 
dom Poisson distributions of N = 2 x 10 4 Ising spins on 
a plane interacting according to Eq. Q. The thermody- 
namic average, M of the polarization for each realization 
was computed by the cluster Monte Carlo algorithm [l4| . 
Results of this computation presented on Fig. 2 suggest 
that 

T c a 13 J, (5) 

i.e. factor C in Eq. is C « 1. 

It should be noted that, strictly speaking, Eqs. Q and 
([5| only apply in the absence of charge carriers. If doping 
is heavy as, e.g., in graphene on SiC 15 , it suppresses the 



transition. Provided that there are more carriers than ad- 
sorbents the Fermi wave number kp exceeds inverse mean 
distance between neighboring adsorbents kp > ^fp and 
conventional Friedel oscillations 
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with the wave vector 
2kp dominate the RKKY interaction in Eq. Q. Under 
these conditions one should expect no ordering. 

Two important remarks on the ordering transition are 
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Fig. 2: Sublattice ordering in the ensemble of 2 x 10 4 adsor- 
bents simulated numerically using cluster Monte Carlo algo- 
rithm [14] , where sublattice occupancies of adatoms are char- 
acterized by the 'Ising spins' rru = ±1 and their thermody- 
namic average, M is computed for 10 realizations of random 
Poisson distributions of spins on a square with periodic bound- 
ary conditions interacting according to Eq. Q. The critical 
temperature, T c ~ 13J, is determined approximately from the 
sharp change in the 'polarization' M. A small, but finite value 
of M at T > T c is a typical finite-size effect. 



due here. One has to deal with competing interactions 
between the adsorbents. In addition to Ising interaction 
in Eq. Q between the adsorbents, each of them causes 
Friedel-like oscillations of a different type: with the wave 
vector equal to the distance between non-equivalent BZ 
corners K and K'. The amplitude is decaying as inverse 
cube of the distance from the adsorbent. From the sym- 
metry point of view these oscillations can be regarded as 
a charge density wave superlattice with a supercell three 
times as big as the unit cell. They generate an effective 
interaction between adsorbents which favors AB configu- 
ration (rather than AA or BB ones) of the closest neighbor. 
Such interaction is anisotropic: when projected onto one 
supercell, two adsorbents determine one of the three pos- 
sible directions of A-B bonds. Due to this anisotropy for a 
low-density p <C aT 2 random coverage the supercell corre- 
lations are unavoidably frustrated and would not lead to 
ordering. This behavior should be contrasted to that of 
adsorbents which reside in the middle of honeycomb lat- 
tice hexagons (such as alkali atoms) or adatoms residing 
over the C-C bonds: such adsorbents cannot distinguish 
between two sublattices, but they can establish partial or- 
dering into a superlattice with a triple-size Kekule-type 
unit cell 
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Another remark concerns the kinetics of the transition. 
The ordered state arises only if adsorbents can hop along 
a flake. This is possible only if the activation barrier for 
the adsorbent's hop between A and B is less than the des- 
orption barrier. Recent studies of hydrogen on graphite 
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suggest that desorption barrier for H on graphite is 
lower than its diffusion barrier. Therefore hydrogenated 
graphene 19 may not be the likeliest candidate to ob- 
serve the suggested sublattice ordering. The alternatives 
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are represented by halogens, e.g. fluorine. The best 
atom/group for realization of the proposed ordering is yet 
to be identified. 

Let us discuss electronic properties of the graphene 
with partially sublattice ordered adsorbents. Opening 
of the spectral gap offers an attractive route to control 
the transport properties. It is quite likely however that 
the Ising-type spontaneous symmetry breaking would re- 
sult in splitting of the flake into domains of opposite "po- 
larization" M = ±1 — regions where adsorbents occupy 
a particular (A or B) sublattice. Such a domain can 
be viewed as a 'weak' topological insulator 20-22 . In- 
deed, while inside each domain the electron spectrum is 
gapped, the boundary between them supports a mid-gap 
one-dimensional band of chiral electronic states 23 - do- 



main wall states (DWS), which are similar to the quan- 
tum Hall edge states [24]. Even in the absence of mag- 
netic field DWS electrons in one valley, e.g. K propagate 
only in one direction, whereas in the other valley, K' the 
DWS on the same boundary propagate in the opposite di- 
rection. Without intervalley scattering, the DWS would 
make a domain wall an ideal ballistic 'wire', with the width 
wdw ~ hv/A ~ (hv / u)(a 2 p)^ 1 and conductance e 2 /h. 

Note that the sublattice-ordered adsorbents still scatter 
electrons between the valleys, at a rate 



pu 2 a 4 "f(e)/h, 7(e) = 2ir/eh z v 



2„,2 



(6) 



Interaction between the oppositely propagating one- 
dimensional DWS leads to their localization, at the length 
scale of the order of the mean free path I for intervalley 
scattering. For the energies we can substitute e ~ A, Eq. 
([!]) into Eq. ^ and estimate the mean free path of a DWS 
as 

hv (hv/ua) 2 



I 



(a 2 P y 



» w DW - 



An immediate consequence of this would be the high re- 
sistivity of graphene flake, with characteristic for hopping 
resistivity temperature dependence, R cx A exp{(T*/T) x } 
(x = g for the Mott hoping law [25], and x = \ for the 
Efros-Shklovskii law [26J ) . As long as I exceeds a typical 
domain size, this system should also display a strong pos- 
itive magnetoresistivity R{B/B*) over a broad magnetic 
field range, ~ h/el 2 . The above-described modification 
of transport characteristics of graphene accompanied by 
opening of a gap in its optical absorption spectrum would 
be natural manifestations of self-organization of a dilute 
ensemble of on-site adsorbents into a sublattice-ordered 
state. 
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